Ab initio calculation of intrinsic spin Hall effect in semiconductors 
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Abstract 

Relativistic band theoretical calculations reveal that intrinsic spin Hall conductivity in hole- 
doped archetypical semiconductors Ge, GaAs and AlAs is large [~ 100(h/ e)(Qcm)~ 1 ], showing the 
possibility of spin Hall effect beyond the four band Luttinger Hamiltonian. The calculated orbital- 
angular-momentum (orbital) Hall conductivity is one order of magnitude smaller, indicating no 
cancellation between the spin and orbital Hall effects in bulk semiconductors. Furthermore, it is 
found that the spin Hall effect can be strongly manipulated by strains, and that the ac spin Hall 
conductivity in the semiconductors is large in pure as well as doped semiconductors. 

PACS numbers: 71.20.-b, 72.20.-i, 72.25.Dc, 85.75.-d 



Spin-orbit coupling is known to be responsible for a whole range of interesting phenomena 
in magnetic materials [lj], including magnetic anisotropy, optical dichroism, and anomalous 
Hall effect. Spin generation and transport in paramagnetic materials can also be induced 
by an electric field because of spin-orbit coupling even in the absence of a magnetic field 
as demonstrated in a number of recent experiments Q. This offers the exciting possibility 
of pure electric driven spintronics in semiconductors, where spin-orbit coupling is relatively 
strong and which can be more readily integrated with conventional electronics [3j]. Earlier 
theories on electric spin generation and transport were based on extrinsic effects due to spin- 
orbit coupling in scattering processes |4]. Recently, it has been proposed that a transverse 
spin current response to an electric field, known as the intrinsic spin Hall effect, can also occur 
in pure crystalline materials due to intrinsic spin-orbit coupling in the band structure 0, |(| . 

A large number of theoretical papers have been written addressing various issues about 
the intrinsic spin Hall effect. In pj, a systematic semi-classical theory of spin transport 
is presented, resolving a discrepancy between the prediction of [5] and the Kubo formula 
result. In jgj, an orbital- angular-momentum (orbital) Hall current is predicted to exist in 
response to an electric field and is found to cancel exactly the spin Hall current in the spin 
Hall effect. There is also an intensive debate about whether the intrinsic spin Hall effect 
remains valid beyond the ballistic transport regime jjj. On the other hand, experimental 
measurements of large spin Hall effects for the Rashba two dimensional electron gas and for 
n-type bulk semiconductors have just been reported Q, till, although more work is needed 
to firmly establish the intrinsic or extrinsic nature of the results. It is expected that similar 
experiments on p-type semiconductors will also appear soon. 

Inspired by the prospect that the intrinsic spin Hall effect may provide a useful prop- 
erty for spintronics applications which is designable based on material parameters, we have 
carried out ab initio calculations on this effect and related phenomena in the archetypical 
p-type semiconductors Si, Ge, GaAs and AlAs. We focus on the p-type because the spin- 
orbit coupling is much stronger in the valence bands, where the intrinsic effect has a better 
chance of dominating over the extrinsic scattering effects. Such results are clearly needed 
to provide a stage for systematic and quantitative comparison between theoretical predic- 
tions and experimental measurements. In addition, the calculations may also help resolving 
some of the theoretical issues for such semiconductors whose discussions are currently all 
based on the four-band Luttinger Hamiltonian for the holes and often with the spherical 
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approximation. 

Our results cover a large range of hole concentration which is beyond the validity regime 
of the four band model. By including all the relevant bands, we find a pronounced spin Hall 
conductivity for all the semiconductors except the light element Si and a roughly quadratic 
dependence of the spin Hall conductivity on the spin-orbit gap. The vanishing of spin Hall 
conductivity in the limit of zero spin-orbit coupling differs qualitatively with the prediction 
of the four band model, which yields a finite spin Hall conductivity even when spin-orbit 
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coupled terms in the Luttinger Hamiltonian approach zero P, 
orbital Hall effect, and find it an order of magnitude weaker than the spin Hall effect, a 
result due to the orbital quenching by the cubic anisotropy of the crystals. This is in sharp 
contrast with the case of the Rashba two-dimensional electron system and the spherical four- 
band model, where exact cancellation between the spin and orbital Hall effects occurs due to 
rotational symmetry in such models |8]. In addition, we calculated the effect of strain, which 
is commonly present in semiconductor multilayers and superlattices due to lattice mismatch, 
and which may be used to enhance or reduce the spin Hall effect. Furthmore, this result on 
strain effect will help the experimentalists to distinguish intrinsic from extrinsic aspect of the 
spin Hall effect in the future [llj|. Finally, we also calculated the ac spin Hall conductivity 
and find it to be large in both pure and hole-doped semiconductors. 

The intrinsic Hall conductivity of a solid can be evaluated by using the Kubo formula 
approach Q3j. In this approach, the intrinsic Hall effect comes from the static cu = limit 

nn 

of the off-diagonal element of the conductivity tensor |y, [13J : 

a xy (u) = --TT (An - /kn') 

< k.n\j x \kn' >< kn'\v y \kn > ... 

Ckn - Ckn' +flUJ + ir] 

where V c is the unit cell volume, Two is the photon energy, and \kn > is the nth Bloch 
state with crystal momentum k. Since all the intrinsic Hall effects are caused by spin-orbit 
coupling, first-principles calculations must be based on a relativistic band theory. In this 
case, the current operators j are — ev, |{/3E z ,v}, and |{/3L 2 ,v} for the anomalous, spin, 
and orbital Hall effects, respectively. Here /3, £ are the well-known 4x4 Dirac matrices [14], 
and v is the velocity operator projected onto states above the electron-positron mass gap 
jis| . Setting n to zero and using Im[l/(x — in)] = tt5(x), one obtains the imaginary part of 
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the off-diagonal element 

<» = ^EE(i-W (2) 

Im[< kn\j x \kn' >< kn'\v y \kn >]<5(e n / n — UuS) 

where e n > n = e^ n i — e^ n . As in previous magneto-optical calculations [l^, we first calculate 
the imaginary part of the a xy and then obtain the real part from the imaginary part by a 
Kramers-Kronig transformation 

<M = -P/ do/ 2 xy[ (3) 

V 7T Jo Uj' 2 - UJ 2 

where P denotes the principal value of the integral. We notice that the anomalous Hall 
conductivity .,(0) of bcc Fe ca.co.ated this waylQ is in good q nan«a t ive a^ent 
with that calculated directly by a generalized Boltzmann equation approach that accounts 
for the Berry phase correction to the group velocity [lfij ]. 



TABLE I: Experimental lattice constant a (see |20| and references therein), average atomic sphere 



21] and references therein) of the semiconductors studied. The 



radius Rw S and band gap E g (see 
calculated band gaps Et he and spin-orbit splitting A so of the top valence bands at T are also listed. 



a (A) R 


ws (a.u 


) E g (eV) Ef e (eV) A. 


w (meV) 


Si 5.431 


2.526 


1.17 


0.81 


47 


Ge 5.650 


2.632 


0.74 


0.28 


278 


AlAs 5.620 


2.615 


2.23 


1.52 


301 


GaAs 5.654 


2.632 


1.52 


0.76 


336 



The relativistic band structure of the semiconductors is calculated using a fully rela- 
tivistic extension Q| of the well established all-electron linear muffin-tin orbital (LMTO) 
method [18]. The calculations are based on density functional theory with generalized gra- 
dient approximation (GGA) Q|. In the fee diamond or zincblend structures, two atoms sit 
at (0,0,0) and (1/4,1/4,1/4), respectively. In the present calculations, two "empty" atomic 
spheres are introduced at the vacant sites (1/2,1/2,1/2) and (3/4,3/4,3/4), respectively, to 
make the structures more close-packed. The four "atoms" in the unit cell are assumed to 
have an equal atomic sphere radius. The experimental lattice constants for Si, Ge, GaAs 
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FIG. 1: (color online) (a) Relativistic band structure of GaAs: (a) Energy bands and (b) density 
of states (DOS) (states/eV/cell). The zero energy is at the top of the valence bands. The dotted 
lines just below the top of the valence bands denote the Fermi levels for the hole concentration 
= 0.025, 0.05, 0.1, and 0.2 e/cell, respectively. 



and AlAs used are listed in Table I together with the corresponding atomic s phe re radii. 
The basis functions used for all the "atoms" are s, p, and d muffin-tin orbitals 18j . In the 
self-consistent electronic structure calculations, eighty-nine /c-points in the fee irreducible 
wedge (IW) of the Brillouin zone (BZ) were used in the tetrahedron BZ integration. The 
calculated band gaps E g and spin-orbit splitting A so of the top valence bands at V for all 
the semiconductors are listed in Table I. As an example, the calculated band structure of 
GaAs is displayed in Fig. 1. 

In the Hall conductivity calculations, a much finer fc-point mesh is needed. Moreover, 
a larger IW (three times the fee IW for Si and Ge, and six times the fee IW for AlAs 
and GaAs which have no spatial inversion symmetry) is necessary for the Hall conductivity 
calculations. The number of fc-points in the IW used is 49395 for Si and Ge and 98790 
for GaAs and AlAs. These numbers are obtained by dividing the TX line into 56 intervals 
(see Fig. 1). Test calculations using 139083 /c-points (80 divisions of the TX line) for Ge 
indicate that the calculated spin Hall conductivity converge to within 2 %. The imaginary 
part of the Hall conductivity is calculated up to 50 eV to ensure that the real part of the 
conductivity obtained through the Kramers-Kronig transformation converges as well. As is 
well known, all the calculated band gaps are smaller than the corresponding experimental 
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values though the calculated spin-orbit splittings A so agree quite well with the experimental 
ones as well as the full-potential calculations (see Table I and |22j). To remedy this defect 
of the GGA, the so-called scissor operator is applied in the conductivity calculations, i.e., 
all the conduction bands are shifted upwards such that the theoretical band gap is the same 
as the experimental one. 
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FIG. 2: (color online) Calculated (a) spin {o s X y) and (b) orbital (<r°y) Hall conductivity as a 
function of hole concentration (n^) for Si, Ge, GaAs and AlAs. (c) Calculated spin Hall conductivity 
for 1%^ = 0.05 e/cell vs. the top valence band spin-orbit splitting (A so ). The dashed line denotes 
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the parabola fitted to the calculated data. The x is the spin-orbit coupling scaling factor 
x = 1 means full spin-orbit coupling and x = means no spin-orbit coupling, = 0.1 e/cell is 
equivalent to n h of 2.5 xlO 21 , 2.2 xlO 21 , 2.2 xlO 21 and 2.3 xlO 21 cm" 3 for Si, Ge, GaAs and AlAs, 
respectively. 



Fig. 2 shows the calculated spin (cr^y) and orbital (cr°y) Hall conductivities as a function 
of hole concentration. The spin Hall conductivity in Ge, GaAs and AlAs increases sharply 
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as the hole concentration rih increases at small rih (rih < 0.02 e/cell) (Fig- 2(a)), being 
consistent with the recent analytical Luttinger model predictions of n]/ 3 [J Q]- However, 
it becomes more or less flat when rih is further increased (Fig. 2(a)). This is due to the fact 
that the spin Hall effect results predominantly from the spin-orbit split heavy-hole (HH) and 
light-hole (LH) pockets centred at T (Fig. 1) |5j. Fig. 1 shows that when rih > 0.05 e/cell, 
the Fermi level is already below the top of the lower spin-orbit split valence band. However, 
there is no substantial increase in the spin Hall conductivity as rih goes through 0.05 e/cell, 
indicating that the lower spin-orbit split valence band has no significant contribution to the 
spin Hall conductivity. Remarkably, the spin conductivity for Ge, GaAs and AlAs is very 
pronounced, being about 100 h/(eflcm) (Fig. 2). These values are comparable to the charge 
conductivity for lightly hole-doped Ge, GaAs and AlAs, in agreement with the prediction 
of the large intrinsic spin Hall effect in hole-doped semiconductors based on the Luttinger 



Hamiltonian 



241 ] . However, our results differ from the predictions of the four band model 



in at least one important aspect. Fig. 2(c) shows that the calculated spin Hall conductivity 
at, e.g., rih = 0.05 e/cell, approaches zero nearly quadratically as the spin-orbit splitting A so 
and hence the spin-orbit coupling strength goes to zero. In contrast, the four band model 
yields a finite spin Hall conductivity even as the spin-orbit coupled terms in the Luttinger 
Hamiltonian approach zero. 
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FIG. 3: (color online) Calculated spin Hall conductivity (a%.y ) (squares) and energy splitting 
(diamonds) of heavy-hole and light-hole bands at T of GaAs as a function of uniaxial strain at the 
hole concentration = 0.1 (e/cell). 



The calculated orbital Hall conductivity is at least ten times smaller than the spin Hall 
conductivity (Fig. 2), demonstrating that the orbital Hall effect in bulk semiconductors is 
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strongly quenched by the cubic crystal fields and covalent bonding. Furthermore, the orbital 
and spin Hall conductivities generally have the same sign when the hole concentration is 
small, i.e., < 0.1 e/cell (Fig. 2). Therefore, the exact cancellation of the intrinsic spin 



Hall effect by the intrinsic orbital Ha 
Rashba spin-orbit term predicted in 



1 effect in the two-dimensional electron gas with the 
would not occur in bulk semiconductors. 
In device applications, semiconductors are typically grown epitaxially on substrates with 
similar lattice constants, resulting in semiconductor multilayers and superlattices. The semi- 
conductor layers are thus generally strained due to small lattice mismatches. The lattice 
mismatch strain may have significant effects on the electronic structure of the semiconduc- 
tors, and hence also on the spin Hall effect j^. In particular, the degeneracy of the HH and 
LH bands at T would be lifted. We have calculated the band structure and Hall conductivity 
of GaAs under [001] uniaxial elastic strains. The calculated spin Hall conductivity as well 
as the energy splitting A of the HH and LH bands at T are displayed as a function of strain 
in Fig. 3. Here the strain e is defined as a = oq(1 + e) where a and ao are the strained 
and unstrained lattice constants perpendicular to [001], respectively. As expected, the A is 
roughly proportional to the strain size. Remarkably, the spin Hall conductivity varies almost 
linearly with the strain, increasing with tensile change and decreasing with compression in 
the lateral directions (Fig. 3). This linear and sensitive dependence on the strain may be 
used to tune the spin Hall effect. 

The shape of the calculated ac spin Hall conductivity spectra for all the semiconductors 
studied here look rather similar, though the peak energy position and peak magnitudes 
may differ. Thus, we only display in Fig. 4 the ac Hall conductivities for GaAs as a 
representative. Interestingly, both the real and imaginary parts of the spin Hall conductivity 
in both pure and hole-doped semiconductors are large (Fig. 4a-b). This suggests that one 
could generate ac spin current in semiconductors without using magnetic field or magnetic 
materials. In contrast, the orbital Hall conductivity in the pure semiconductors is zero 
(i.e., within numerical noise) whilst it is two order of magnitude smaller than the spin Hall 
conductivity in the hole-doped semiconductors (Fig. 4). 
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and National Science Council of ROC, NSFC (No. 10404035) of PRC, US Department 
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FIG. 4: (color online) Calculated spin and orbital (o~°y) Hall conductivity of GaAs as a 

function of frequency at the hole concentrations = 0.0,0.1 (e/cell), respectively. 
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